Uniform asymptotics of the coefficients of 
unitary moment polynomials 

By G.A. Hiary and M.O. Rubinstein f 

Keating and Snaith showed that the 2k*^ absolute moment of the characteristic 
polynomial of a random unitary matrix evaluated on the unit circle is given by a 
polynomial of degree k^. In this article, uniform asymptotics for the coefficients 
of that polynomial are derived, and a maximal coefficient is located. Some of the 
asymptotics are given in explicit form. Numerical data to support these calculations 
are presented. Some apparent connections between random matrix theory and the 
^~^ Riemann zeta function are discussed. 
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(N 
5P 1. Introduction 

■^ Let Zn{A,9) denote the characteristic polynomial of an N x N unitary matrix 

A evaluated at exp{i9). Let 9n denote the eigenphases of A. Then, Zn{A,9) := 
det (e*^/ — A) = rinCs*^ — e'^"), Keating and Snaith [KSj gave an explicit expression 

I— I for the expected moments of \Zn{A, 9)\: 

^ KN\ZAA,9)r^f[W^^^, 5Rfc>-l/2. (1.1) 

d 

C The expectation is taken with respect to the normalized Haar measure on the 

1—^ group oi N X N unitary matrices U{N). This can also be written in terms of 

the Barnes G-function, which is the entire function of order 2 defined by G{z + 
1) = (27r)^/2e-(^+(i+7)^')/2 ]-[«=^^(i_^2/„)«g-^+zV(2")^ and satisfying the relations 
G(l) = 1, G{z + 1) = T{z)G{z). The r.h.s. of ([TT]) equals 

O 

^ G{k + If G{N + 1)G(A^ + 2fc + 1) 

^^ G(2fc + 1) G(Ar + fc + l)2 ■ ^ ' ' 

C^ When 2k G Z, this simplifies to a polynomial in N of degree k^. This polynomial 

f^ can be expressed in several equivalent ways, each of which is useful: 

> E^|Z^(A,0)|2'=^gg_Lll! Yl iN + ^+J + l) (1.3) 

'"^ ^ ' 0<i.j<k-l 



cn 



oo 
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-'■-'- (7 + fc)! J--'-' ■'^ 11 ^ •" 11 V i + j + 1 

j=0 ^-^ ' j = l j=k+l Q<i,j<k-1 ^ ■' 
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Expanding, we define, 



(1.4) 



r=0 



The asymptotics of tfie leading coefficient co(fc) can be obtained from tlie asymp- 
totics of tlie Barnes function: logG(z + 1) ~ z^(log(z)/2 - 3/4) + log(27r)z/2 - 
log(^)/12 + C'(-l) + 0(l/2), from which it follows that co(fc) = -fc^ logfc-fc^ log4+ 
3fcV2 - log(fc)/12 + C'(-l) + 0{l/k). 

(a) Results 

The main purpose of this paper is to develop uniform asymptotics in r for the 
coefficients Cr(fc), < r < fc^, as A: — )■ oo. We also obtain some explicit asymptotics, 
and estimate in r, for given fc, the maximal Cr{k). 

We prove the following theorem in the next two sections as a sequence of five 
lemmas: 



Theorem 1. Let u := Ur^k be the unique positive number satisfying u'^^^^ - 



+'j 



'"Ej=fc+i -^ == P ~ r. Further, define U := Ur^k by U = u^^j^i (^^qrj)2 



X)j=fe+i ''(u+ Xi ■ Then, uniformly in r ^ (0, k ), and as k -^ oo, 



^(^)-©' 



fc'\ Pk{ 



u 






fc^-r+1/2 



.1/2 



1 + 



logfc 



Furthermore, if we fix a < 1, then, uniformly in < r < k"' 



co{k) r\ 



l + 0(fc2(°-i)) 



l + 0(/c2("-i)) 



(1.5) 



(1.6) 



and 



Cfc2_r(fc) 



[kAy 



l + 0(fc2("-i) logfc) 



fc2\ A"- 



l + 0(fc2("-i) logfc)] (1.7) 



with the implied constants in the O terms depending on a and where 

2k 

A = 2 ^ l/j-^log4+- + 0(l/fc2) 



j=k+i 



(1.8) 



Regarding the location of the maximal Cr{k) we prove in section [5] 
Theorem 2. Let 



2 k 



Ai 



2k -J 



E-XT+ E --T^ = fclog4-log(fc/2) + l/2-7 + 0(l/fc), (1.9) 



j=i 



j=k+i 



where 7 = 0.57721 . . . is the Euler constant. Then, there exists p > such that, for 
all k sufficiently large, a maximal Cr{k) occurs for some 



re[k^-fi- p\og{kf/k, fc2 - /i + 1 + p\og{kf/k] 
and no maximal Cj.(fc) occurs outside of that interval. 



(1.10) 
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Notice that the size of the interval in (1.10 1 is shghtly larger than 1. So typically, 
it will only contain one integer, and the maximum with occur for r = \k'^ — fi\ , but 
it can contain two integers in the event that /x is very close to an integer. 

In section [4] we derive a more precise formula for the leading coefficients, valid 
for r = 0{k°') with a < 2, the first few terms of which are listed below: 

r(r-l)(583r2-715r + 183) 
6480fc6 



Theore 


m 3. 


Cr{k) 

co{k) 


P"" / 7r(r-l) r(r-l)(26r-15) 
r! ^^^^l^ 12fc2 144/fc4 




r{r ~ l)(2758r3 - 4499r2 + 463r + 1491) 



51840A:8 ' ' ' fc2j 



(1.11) 



The j-th term in the exponent is of the form qj{r)/k'^^ , where qj{r) is a polynomial 
in Q[r] of degree j + 1, divisible by r{r — 1), and satisfying the bound qj{r) = 
0{jlog{j)r{\r)^), for some A > 0. 

The above expansion can be transformed into a similar one involving the bino- 
mial coefficient ( ^ ) : 

Crik) _ ,rfk^\ ( r{r-l) r(r-l)(2r-3) r(r - l)(43r2 - 175r + 183) 
'^^ ^ V''j'^'''^V 12^?^ 144fc4 6480fc6 

r(r- I)(166r3-611r2 + 31r + 1059) qj(r) \ 

51840fc8 ' ■ ■ " fc2j -•••]• 



(1.12) 



This can be achieved by writing C; ) ^ = J]^^! (l - ^x) = exp(- Yl'jJi Em=i :^£^) 



exp(- Y.m=i ^'"tm+i)mk^"^ )^ whcrc Bjnix) are the Bernoulh polynomials, and we 



have used (1.18). 

The apparent smaller coefficients in the exponent of ( 1.12 ) in comparison to (|1.11[) 



suggests that the binomial coefficient appears naturally as a factor in the asymp- 
totics. A heuristic explanation is provided in the next section. 

It turns out the problem of deriving uniform asymptotics for the Cr(A;)'s is similar 
to that of deriving uniform asymptotics for the unsigned Stirling numbers of the 
first kind. These numbers are usually denoted by S'(fc, r), and can be defined by the 
generating function Oi^o (^ + .?) ~- ^r=i '^(^^ r)x^ . The above generating function 



bears some resemblance to the generating function of the Cr{k)^s given in (1.3) 
which is the reason the settings in the two problems are similar. 

Our proof of lemma Is] parallels the proof of asymptotic (1.3) in Moser and 
Wyman |MW| . where uniform asymptotics for the S{k,r)'s were developed. See 
also Theorem 2 in the later work of |CRT| . which relies on the work of [MWj . 

(&) Motivation and connections to the Riemann zeta function 

For a fixed k with 3fffc > —1/2, and as T — )■ 00, the 2A:*'' moment of the Riemann 
zeta function on the critical line C(l/2 + it) is conjectured to have an asymptotic 
expansion: /^ |C(l/2 + it)\'^'^ dt ^ /^ P/i;(log(i/(27r)))(it, where Pk{x) is an asymp- 
totic series in descending powers of x, of the form: Pk{x) :— '^'^QCr{k)x ^^. 
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When X e Z, then this series terminates and is a polynomial of degree fc^ in x: 

The leading coefficient is expected to be of the form co(A;) — akgk/k'^^-, where 
ttfe is a certain, generally understood, "arithmetic factor," and gk is an integer 
for integer values of k (see |CGlj ). It is well-known gi = 1, g2 = 2, and it was 
conjectured g^ = 42, 34 = 24024, in |CG2] and |CGo| . respectively. But little was 
known about gk in general. 

Based on the analogous computation in random matrix theory, Keating and 
Snaith [KS] conjectured that gk/{k^l) = G'^{k + 1)/G(2A: + 1), or co(fc) = co(fc)/afc. 
In other words, the leading coefficient oi Pk{x) and of Pk{x) are expected to coincide, 
except for "arithmetic effects." 

Recently, Conrey, Farmer, Keating, Rubinstein, and Snaith [CFKRSl] conjec- 
tured, for integer fc, a formula for the polynomial P expressed in terms of 2fc-fold 
residue. The residue formula allows one to derive complicated formulae for the 
Crikys |(]FKES2j . 

Precise information concerning the extreme values of the zeta function up to 
given height could be derived from precise knowledge of the uniform asymptotics 
of its moments in both the k and T aspects |FGH| . However, it is not clear how to 
obtain the uniform asymptotics from the formulae in [CFKRST] and |CFKRS2] . 

As a first step, we decided to try and understand the behaviour of the coefficients 
of the moment polynomials. We hoped that deriving asymptotics for the coefficients 
that occur on the random matrix theory side (i.e. the Cr(fc)'s) could shed some light 
on the zeta function side (i.e. the Cr(fc)'s). In a sense, that turned out to be the 
case. In an upcoming paper |IIR| . we show for a fixed a < 1, and uniformly in 
< r < fc", it holds, 

coik) 

where 7 = 0.5772 ... is Euler's constant, and Bk is a certain explicit number 
(see [HRj for details) satisfying Bk ~ 2fclogA:. To facilitate comparison, we give 
the corresponding result for unitary random matrices (see lemma pi). 



{2Bk + 2"fkyrj l + 0(/c2("-i)) , asfc^oo, (1.13) 



Co{k) 



l + 0(A:2("-i)) , asfc^oo. (1.14) 



(c) Basic facts 

For a fixed r, one can easily show that Cr{k)/co{k) is a polynomial in k of degree 
3r with leading coefficient 1/r!: 

Furthermore, pulling out the 1/r!, all the coefficients of rlbr{k) are polynomials 
themselves in Q[r]. These facts can be proven by writing Cr{k), which are express- 



ible from (1.3) as elementary symmetric polynomials, in terms of power sums via 
Newton's identities applied to the set of k^ numbers {1, 2, 2,3,3,3, ... ,2k — 2, 2k ~ 
2,2/c-l}: 

r 
rbrik) = Y.(-ir~^p^{k)br-n{k), (1.16) 
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where Pn{k) are the sums of powers: 

k 2k 

Pn(fc)=Ej^'"+ E (2^-j)j"- (1-17) 

j=l j=k+l 

The above can be expressed in terms of Bernoulli polynomials using 

fc 






'^-^ m + 1 



Thus, breaking (1.17) into three sums, we find 



.,. _ Bn+2ik + 1) - Bn+2{1) , ^,^ g„+l(2fc)-B„+i(fc+l) g„+2(2fc)-g„+2(fc + l) 

^"^^^ " ^T2 ^^'^ ^n ^+2 

2E„+2(fc+l)-5„+2(2fc)-g„+2(l) ^ ^^^ gn+l(2fc)-g„+i(fc+l) 

n+2 n+ 1 

on+2 _ 9 

^ ^ ■.n+2 



fc"+^ + ... (1.19) 



(n + l)(n + 2) 

is a polynomial in k of degree n+2 and leading coefficient (2"+^ — 2)/((ri+l)(n+2)). 
The last step follows from the expansion of Bn{x) in terms of Bernoulli numbers. 



Therefore, from the n — 1 term in (1.16), we inductively get the leading term 



of (1.15), and the recursion also gives the coefRcients of br{k) as rational numbers 
(functions of r). We fist the first eight p„(fc): pi(fc) — fc"^, P2{k) = 7/6 k^ — 1/6 A:^, 
P3(fc) = 3/2 k^ - 1/2 k^, pi{k) = 31/15 fc6 - 7/6 fc^ + 1/10 k^, pz{k) = 3 F - 5/2 k^ + 
1/2 fc3, pe(fc) = 127/28 fc^ - 31/6 fc^ + 7/4 k^ - 5/42 fc^, p7(fc) = 85/12 fc^ - 21/2 fc^ + 
21/4 fcS _ 5/6 fc3, pfi{k) = 511/45 /fcio - 127/6 fc« + 217/15 fc^ - 35/9 fc* + 7/30 fc^, 
and the first few ^^(fc): 6i(fc) = fc^, 52(fc) = 1/2 fc*^ - 7/12 fc* + 1/12 fc^, b3,{k) = 
1/6 fc^ - 7/12 fc^ + 7/12 fcS - 1/6 fc3, hi{k) = 1/24 fc^^ _ 7/24 fc^o + 205/288 fc^ - 
527/720 fcS + 85/288 fc* - 1/40 fc^, 65(fc) = 1/120 fc^^ - 7/72 fc^^ + 125/288 fc" - 
677/720 fc9 + 1489/1440 fc^ - 97/180 fc^ + 1/10 fc^, be{k) = 1/720 fc^* - 7/288 fc^^ + 
11/64 fci* - 32927/51840 fc^^ + 22931/17280 k^° - 38245/24192 fc« + 10513/10368 fc^ - 
47/160fc* + 5/252 fc2. 

One readily observes from the powers of fc that appear that br{—k) = {—lYbr{k), 
and this can be proven inductively from the recursion (first establishing this prop- 
erty for Pn{k) from its expression in terms of Bernoulli polynomials). 

The recursion also allows us to work out specific formulas for the lower terms 
of br{k). For example, the next to leading term has degree 3r — 2. Writing 

br{k) = Y^ br^jk^''-^^ (1.20) 

0<j<(3r-2)/2 

we have, from the n = 1,2 terms of the recursion, the relation: rbr^i = &r-i,i — 
7/6br-2,o- But 6r-2,o = l/(^ ^ 2)!, and one easily checks inductively that 6^,1 = 



12 (i--2)! ■ 

More generally, writing p„(fc) = X)o<)<= Pn-j^"^^^^"'' ^^ have, on comparing 



the coefficient of fc'''' ^^ on both sides of (1.16), that 



j+l j+l-n 

^br,j = y^(~l)"~^ y^ Pn.abr^n,j + l-n-a- (1-21) 

n—1 a— 
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One can show, inductively on j, that h^j is of the form 



br,3 = -^ (1-22) 



where gj (r) is a polynomial in r with rational coefficients, by plugging this into ( 1.21 ), 
and separating the (n, a) = (1,0) term from the rest as follows. Putting the 
other terms over a common denominator, and using pi o — 1, we have: r ^^ , — 
^(r-i)\ + (r-iv. ' ■where, hj (r) is the polynomial obtained by our inductive hypoth- 
esis from the terms {n,a) ^ (1,0). Therefore gj{r) — gj{r — 1) = hj{r), and this 
difference equation allows us to solve for all but the constant coefficient of gj , which 
can be obtained from a specific value of brj ■ 

For example, to work out 6^2, we substitute 6^.2 = 52(^)/^! into the recursion 
rbr,2 = Y^n=ii-^T~^ Y^a=o Pn,abr-n,3-n-a, and puU out the {u, a) = (1,0) term, 
giving, r^ = ^fey!^ + 72(r-4)! + 6(?=2)I + 2ii^. ' Clearing denominators we 
have 52(0 - g2{r - 1) = ||(r - l)(r - 2)(r - 3) + i(r - 1) + |(r - l)(r - 2) = 
49/72 r^ - 31/12 r^ + 227/72 r - 5/4, and solving for the coefficients of 32 yields 
ff2(r) = 49/288r4-25/48r3 + 131/288r2-5/48r+52(0). From &2,2 = 52(2)/2! = 1/12 
we get 52(0) = 0. Hence, after factoring the above, 6^,2 — '^^288(r°-2)! ^° • ^'^ ^^^^ ^^^ 
first few terms in the expansion of br{k): 

_ 1 3. 7 1 3.-2 49r^-101r + 30 3.-4 

'''■^''^ "^'^ "12(7^2)!'' + 288(r-2)! '^ 

1715 r^ - 5460 r^ + 4069 r - 732 3^,6 
51840(r-3)! 

12005 r'^ - 52430 r^ + 69967 r^ - 22726 r - 11928 3^_8 
"^ 2488320(r-4)! 

(117649 r'' - 1130871 r'^ + 3998449 r'' - 6072801 r^ + 2402926 r^ + 383976 r + 6420672) ^^_^ 



209018880(r--4)! 
49076440 r® + 225 
+297681419 r^ - 108063712 r^ + 978267588 r + 1471379040)fc^ 



+ 75246796800(r-5)! ^^^^^^^^ ' " 49076440 r« + 225080030 r' - 466558120 r* 



(1.23) 



2. Basic asymptotics 

The coefficients of Pk{x) can be expressed as elementary symmetric polynomials 
in its roots. For example, if we let —zi, —Z2, . . . , —2^2 denote the roots of P}~{x), 



then, from the last equation in (1.3), we have C)^'i{k) = 1, C}^2_i{k) = X^KiXfe^ 'h^ 
Cfe2_2(fc) = y^Ki ^o <t2 , ^ , and so on. The question is reduced to estimating 

' ^ -'-_Jl '^J2_i"' -^j 1 ^J2 

y^Ki ^1 <- ^> <t2 ~ • If ^ is small, then the terms Zj. in this expression 

may be thought of hcuristically as independent random variables drawn from the 
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distribution 



[z = J) = 



j/k\ iije{i,...,k} 

{2k-j)/k\ iije{k + i,...,2k} 



(2.1) 



Therefore, for small r, one may expect something like the following to hold, Cfc2_^(fc) ; 

C;) [Ep {l)Y = (^ 1^, where A is defined in ([L8|. This is precisely the state- 
ment of our lemma [l] 

Lemma 1. Fix a < 1. Assume < r < A:". Then, uniformly in r over that range, 
and as fc — >■ cx), it holds 



Ck^-r{k) 



{kAy 



l+0(fc2("-l)logfc) 



fc2\ A'- 



l + 0(fc2("-i)logfc)j . (2.2) 



Proof. One may write 



where 



(fc 2k 

logHii+^/jF n i^+^/jf^' 
3 = 1 J = fc+1 



k 2k 

i=i j=fe+i 



(2.3) 



(2.4) 



By the Euler-Maclaurin summation formula, 771 = /clog 4 + 0(1), and 77™ < 21ogfc 
for TO > 1. So, we anticipate the main contribution will come from the term 
exp(77ia;). We thus write 



exp(77ia:;) exp 2_, Vrax"^ =: exp(?7ia;) f{x) , 



\m=2 



and expand / in a Taylor series 

00 
/(x) = l + 5]/3„x" 



A direct application of Cauchy's estimate yields the upper bound 

2m Jc 2"+^ 



\Pn 



dz 



0(l0g(fc)"/2) , 



(2.5) 



(2.6) 



(2.7) 



where C is the circle centered the origin of radius 1/^/Togk, and k sufficiently large. 
Now, Cfc2_r(fc) is the the coefficient of x'' in the Taylor expansion about zero of the 



function (2.5 1. As (3i —0 in expansion (2.6), it follows 

7^ 7" r — -n />-. / T 



Cfe2_r(KJ 



E 



Vi 



r\ ^-^ (r — n)\ 



Pn 



m_ 



1 + E 



r\ P„ 



,r — n)\ n? 

n=2 "• ' '1 , 



(2.8) 
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Fix a < 1, and let re [0, fe°] . Then estimate (2.7) yields Cfc2_^ (fc) = ^ [l + O (fc2("-i) logfc)] , 
with the implied constant depending on a. Since ryi = fcA, we arrive at 

Cfc._,(fc) = ^ [l + 0(fc2("-i) log k)\ , (2.9) 

r! L J 

as required. Finally, to obtain the result in terms of the binomial coefficient, note 



A"- /fc2\ 
k^'Xr J ~ 


[kAY ei 
r! (/c2-r)!fc2'- 


{kAY 
r! 



r-l 

exp|log[](l-j7fc2) 



(2.10) 



and 



exp logn(l-j7fc') =1 + EE 



j=i 



I j=i ™=i 



m fc2 



l + 0(fc2("-i)). (2.11) 



D 



A similar heuristic applied to the leading coefficients of Pk {x) leads one to expect 
that Cr{k) sa co(fc)( ^ ) []Ep(z)]'^ = co(fc)( ^ )A:'', in agreement with lemma 2 below. 

Lemma 2. Fia; a < 1. Assume < r < fc". Then, uniformly in r over that range, 
and as k ~> CO, it holds 



Cr{k) _ fc3'- 

Co{k) r! 



l + 0(fc2("-l)) 



fc'' 



l + 0(/c2("-l)) 



(2.12) 



Proof. The leading coefficients of P^. (x) are the trailing coefficients of a;*^ Pk{l/x) - 
co(fc) ni=i(l + Ja^)-' nj=fe+i(l +.?a;)^''"-'- Also, by Taylor expansions, logJl^^ill 



r2fe 



J^P n;:fe+i(l + J^)''-^ = E™=i ^^^;^P™(fc)x", where p„(fc) = E;=i J 



771 + 1 



^2/c 



Sj=fc+i(2^^i)j™- Now, Cr{k)/ca{k) is the coefficient of x'' in the expansion about 
zero of 



exp(pi(/c)a;)exp ^ 



-1) 



771+1 



m 



-Pmik)x" 



eMPiik)x) {l+J2 ''m{k)x"' (2.13) 



Since Pm(k) < (2A:)™+2 for m > 1, Cauchy's estimate supplies the bound Vm{k) — 
0{k'^"^), as fc — )• cx). Also, it easy to verify pi(fc) = k^ and I'lik) = 0. 

Finally, suppose r € [0, A;"], where a < 1 is a fixed constant. Then, put together, 
we have 



Co(fc) r! ^—^ [r -- m)\ r\ 



(2.14) 



D 



By taking more terms in the above lemma, we can obtain a more precise asymp- 
totic, though still valid only up to r = 0(fc") with a < 1. For example, truncating 
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the sum in (|2.14) at m — 5, gives 



c.(fc)^fc3'-/ (r-l)r(7fc^-l) (r-2)(r-l)r(3fc^-l) 

co(fc) r! i^ 12fc4 + 6A:6 ^ ' 

(r - 3) (r - 2) (r - 1) r (245 k^ - 814 fc^ + 425 k^ - 36) 

1440fcio 
(r - 4) (r - 3) (r - 2) (r - 1) r (105 fc^ - 266 k^ + 185 fc^ - 36) 



O 



(fc6(-i))y 



360fci2 
One can also rearrange this expansion, collecting terms according to the power of 



1/k to get another derivation of the expansion (1.231. 

The situation away from the tails is more complicated. There, we apply a saddle- 
point technique in a similar way to jMWj , where it was used to develop asymptotics 
for Stirling numbers of the first kind. Due to such similarities, some of the details 
in the proof of lemma [3] are only sketched. 

Lemma 3. Fix < a < 2. Let u :— Ur,k be the (unique) positive number satisfying 

uY^ + uY^^^^^k'^r. (2.16) 

Further, define U :— Ur,k by 

TT \^ ^^ , \^ JC^k-j) 

3 = \ ^ •" j = fc+l ^ ' 

Then, for \ < r <k'^ — \, and as k ^ oo, we have 

^'■(^) = j^% r [1 + ^(1/^)] ■ (2.18) 

// we restrict r to lie in the interval fc" < r < k^ — k"^ , then, as k ^ go, 

Notice that because we can take a < 2, the interval covered by this lemma, 
fc" < r < fc2 — /c" , overlaps with the two intervals in lemmas [l] and [2J 

Proof. By Cauchy's theorem, Cr{k) — ^ J^ k^ J+i ^^^ where C is any contour 
circling the origin once in a positive direction, and 

fe 2fc 

Pk{z) = co{k)l[{z + jy lliz+jf''-^. (2.20) 

j=i j=fe+i 

The plan is to obtain very good approximations of the above integral via a saddle- 
point method. So, consider the function, 

f{z) := logPfe(z) - (fc2 - r) log z . (2.21) 
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The saddle points of f{z) are the solutions of 

k . 2fe „, . ,2 

/'(-)-E^+ E ^^-^L^^O. (2.22) 

By monotonicity of z^ ■ -^ j/(j + z) + z^ ■ j,^j^(2A; — j)/(j/ + 2) for z > 0, equation 



(2.22) has a unique positive solution, which we denote by u :— u^k. From here on, 
we follow the standard saddle point recipe (see |DBj ). in a generally similar way 
to what was done in |MW| . The contour C should cross the saddle-point u in the 
direction of steepest descent. The direction of steepest-descent is determined by the 
argument of /"(u). Since, 

/"W = ^->:7-i4^- >: 7^^^-^> (2.23) 



fc2-r 


k 


2k „, 


Z2 



then arg f"{u) — (because uf"{u) > —f'{u) = 0, as can be seen by showing 
that zf'{z) is increasing). So, as in the standard saddle-point recipe, the angle of 
passage should be 7r/2 — (1/2) arg f"{u) = 7r/2. This suggests the contour choice 
C := {ue^^ : — tt < 9 < vr}, which is the same as in [MW| . Writing 

cr{k)^^ r ef^^-'Ue, (2.24) 

27r J_^ 

we have, by calculations analogous to those in |MW| . 

Cr{k) = — / e^^^"^"' d6 + (ef(^)-^"u\ (2.25) 

27rJ|e|<e ^ ^ 

for some absolute constant c > 0, and e > a small number to be chosen later, and 



E(^rT7F + ".5,7^rT7F- ^''''^ 



de 






Next, expand f{ue'^^) about 61 = to obtain (for 6 < 1/100 say) 

oo 

f{ue^') = f{u) - 2 ^' + E >"(")^'" • (2.27) 

In that expansion, we used the fact u is a saddle point, i.e. ^/(ue'^)L — 0. 

Now, Cauchy's estimate applied to 272/(^6*^) yields 7,n(u) = 0{U), for m > 3. So, 
choose e such that, 

^^ljS-1/2^ 0<(5<l/6, J is fixed. (2.28) 

Consequently, 



f(u) r / °° \ 

c.(fc) = ^ y^ ^ e-f ^M 1 + ^ M„0™ \de + (e/(")-^^^) . (2.29) 



Article submitted to Royal Society 



Uniform asymptotics of the coefficients of unitary moment polynomials 11 



Cauchy's estimate, in combination with the previous estimate on the 7m 's, gives 
/im = 0{U"^^^), uniformly in m and U. (Notice by the assumption 1 < r < fc^ — 1 
in the statement of the lemma, combined with lemmaH] we have [/ > 1.5 for all 1 < 
r < P_i and k large en ough , and so |At™6''"| = 0([/-(i/6-<5)m-) ^ 0(i.5-(i/6-5)m)^ 

Therefore, the series in (2.291 converges geometrically even at extreme values of r; 
i.e. for r > 1 and r < k"' — 1.) We can also get a sharper estimate for m fixed. 
Exponentiating, e^"=3''''"^"^^", ensures ^„i = Om(C/L™/'^J). 



So, if the sum over m in expression (2.291 is truncated at some integer M > 3, 
then, using the first bound on /i™. 



Cr{k) 



of{u) 



( ^^-1 \ 



2^ ■he\<^' V ' ^3 



+0 [[/'5-l/2-(l/6-5)Mg/(«)\ _|_ Q LfM-c. 

The domain of integration is extended to (— oo, oo), to obtain 



(2.30) 



.{k)^ 



2n 



e-¥«^ 1 



A/-1 

m— 3 



+0 (e^(" 



)-c(72 



o 



— / e-*'' 1+ V Ai,„rM d0 (2.31) 



The function ^^'"e ^^ /^ achieves its maximum at 6' = y^2m/U. It follows, on 
estimating the integrand, that 



52mg-i 



de = o 



e 2 



U{l-25)m \ ' 



for 2to < [/^^ 



(2.32) 



Thus, if the condition 2AI < U is imposed in equation (2.31 1, we obtain 

/(«) /-oo / Af-1 \ 



27r 



771 — 3 



V2TrU 



L(M-1)/2J 

E 

m=2 



2'"r(m + 1/2) /i2„ 
^/7^C/'" 



+ O fuS-{i/e-S)M\ 



(2.33) 



The first part of the lemma follows by taking S = 1/12, M = 13 say, and the 
estimate fj,m = Om(C/L™/'^J). For the second part we observe that 



U = f7(fc") for k°' < r < k'^ ~ k°' 



(2.34) 



as fc — >■ CX3. This bound can be proven as follows. First, recall, by definition, 
is not hard to see, for x > fc/2, x^-^ , ■+xy ~ ^{k^/^)- Also, for < a; < k/2, 

■ 2 

^'^x<i<k I +x)'^ ~ n{xk). By lemma 4, and the monotonicity of the saddle point 
Ur as a function of r, we have A:""^ << u,. << fc"^^" for k" < r < A;^^". Combined 
with the above two lower bounds, this yields (2.34). D 
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3. A uniform asymptotic 

One unsatisfying feature of lemma [3] is the implicit nature of the asymptotic ap- 
proximation it provides. Also, it might be desirable to unify lemmas IT] through [3] 
into a single asymptotic approximation applicable over all regions. The purpose of 
this section is to resolve these issues. Lemma [4] provides explicit approximations to 
Ur^k and Ur,k- Lemma ^ provides a uniform asymptotic for Cr{k). 

Lemma 4. Fix a < 2. Let Ur.k and Ur,k be as in the statement of lemma[^ and 
let A be defined as in { 1.8), so A = log4+ l/(2k) + 0{\/k^). Then uniformly in r, 



kA 



[l + 0(fc"-2iogfc)] ^ ifre{k'^-k°',k^), 



Urh 



k^/r [1 + 0(/c"-2)j ^ ifre{0,k°'). 



Ur,k — 



(fc2-r) [l + 0(fc"-2logfc)] , if r £ {k"^ - k°' ,k^) , 

r [l + 0(fc""^)] , if re (0,fc"). 

Proof. Recall the saddle point u :— Ur,k is the positive number satisfying 

uY^+uy'^^^^^e-r. (3.1) 



Assume r G (fc^ — fc", fc^). The left-hand side of equation (3.1 ) is strictly increasing 
for u e (0,oo), so Ur^k is unique, and it is easy to see, by direct substitution for 
Ur^k, that 

M^fc = 0((fc2_r)/fc) =0(fc"-i). (3.2) 



Let us rewrite equation (3.1 ) as u^i<^.<2„ t+-j^'^^2u<j<k 1+1^7+" Si=fc+i T^T 



Then, uEi<j<2«j/(" + j) = 0(m2), uE2u<j<fc 1/(1 + "/j) ^ ku + 
0{u'^ logfc), and uEj=fe+i(2'fc/i - 1)/(1 + u/j) = kuA - ku + 0(u^). Therefore, 
ukA + 0{v? log k) = k^ — r, and we see, in combination with ( |3.2[ ) that 

As for the case r e (0, fc"), note if u is replaced by fc^/(10r), the l.h.s of (3.1) is 

^EyTTTrf^ ^ + ^ T , , fnn^ . <k^-lOr + Oir^k') = k" - 10r[l + 0{r/k^)] .(3.4) 

lOr -^ A:3/(10r) +j lOr .^^ fc3/(10r)+j ~ \ i ) l \ i ji \ j 



Therefore, by the monotonicity of the left side of equation (3.1), the saddle point 
satis fies u = ^{k^/r) ~ il.{k^^°'). Using this estimate on u, we can rewrite equation 



(|3J| in the form J2''j=i T^ + T,%k+i TTjJTi =k^-r, which implies ^2"]=! i(l - 
j/'a) + J2j=k+ii'^k~j){l—j/u) + 0{k'^/u^) = k'^ — r. From this, it is straightforward 
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to verify k^/u+O {k^ /v?) = r, or u — k^/r [l + 0(fc" ^)] , as claimed. Asymptotics 
for U :— Ur.k are derived similarly. We have, 



ikA + 0{u^ log k) = {P -r)[l + (/j"-2 log ^)j ^ if r G (fc2 - k°',P) , 



U 



k^/u + O (k^/u^) ^r[l+ (fc""2)] , 



if re (0,fc") 



D 



With the aid of lemma |4J it becomes possible to combine lemmas [T] through [3] 
into a single lemma. 

Lemma 5. Let u :~ Ur.k and U := Ur.k be defined as in the statement of lemma[^ 
Then, uniformly in r (£ (0, k'^), and as fc —)■ oo, it holds 



^ ' fc2! \r ) ^yk^-r^ ' 



\k^J \r J ^u'"'-'- V k^J 



l + O 
1 + 



/logA: 

^fc2/3 

log A; 



f,2/3 



(3.5) 



Note: we exclude r — 0,k from the statement of the lemma because the r.h.s. 
vanishes at those points. 

Proof. Lemmas fll [2] and p^ give the asymptotics for Cr{k) in three overlapping 
regions that cover all < r < fc^. We compare these three asymptotics to the r.h.s. 
above, and find the a in those lemmas that gives the optimal remainder term (it 
will transpire that a ~ 3/4). 

So fix a < 1, and suppose r € (0, fc"). By lemma 111 we have 



u = fc7r[l + 0(fc"^2)], C/ = r[l + 0(fc"-2) 



(3.6) 



=^ logu = log(fcV?-) + 0(fc""^), logC/==logr + 0(A:"-2). (3.7) 

Using log(l + x) = X + 0{x^) with x = j/u, and the above estimate for u, we get 

k 2fc 

log{Pk{u)/co{k)) = ^jlog(j+M)+ ^ {2k-j)log{j+u) 

3 = 1 J = k+1 

k 2k 

= k'logu+^Y.^' + - E (2fc-j)j + 0(fc2"-2) 
u ^-^ u ^-^ 

j=l j=k+l 

2i„„.. , ;,3 /„, , ^/,,2q-2 



= k^ log u + kyu + oik^"-^) . 
And by Stirling's formula. 



r J (fc2 - r)! y/2^ (fc2 _ ^)fe^-r+i/2 



l + O 



fc2-: 



(3.8) 



(3.9) 
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Put together, for r G (0, fc"), 
^e'^^ fk'^\ Pk{u) 



co(fc) 

Co(fc) 



Uu 



k^ — r 



ri^2 _^k^-r+l/2^r+l/2 



,rlog(*:^/r)+r-(l/2) logr+rlogr+(l/2)logr-r+0(fe^° ^) ,,^ K 



3r 



(3.10) 



co(fc)^[l + 0(A:2"-2)] 



Lastly, by lemnia[2] 



Cr{k) 



V27 



,-fe^ 



Pk{u) 



.^k" - r)'='-'^+i/2 ^-+1/2 [1 + o [k"^-^)] , r e (0, A;"X3.11) 



Next, we consider the case r G {k", k^ — fc"). By Sthling's formula 



P! 



(fc^ - ^)fc -'^+l/2j.r+l/2 



~fc2! 27r(/c2-r)'='-'^+i/2r'^+i/2 

[l + 0(fc-")]. 



l + O 



1 1 



fc2 



(3.12) 



So by lemma ^ and for r e (fc", fc^ — fc") 



Cr{k) 



{!^,)^^Sk'--f-'^"'^'^"'[^ + 0{k-'^)] . (3.13) 



Similar manipulations together with the remark following lenimalllgive: for r e (fc^ — fc", fc^), 

C,(fc) = ^^?^f^ f'l ) ^#^(fc^ - rf-r+l/2 ^.+ 1/2 ^1 ^ o (^2.-2 Jog fc)] (3.14) 



P! 



r J VUu'"''-'- 



Finally, comparing (3.11), (3.13), and (3.14), we find that the optimal a satisfies 



2 — 2a = a^a = 2/3, giving the O term claimed in the lemma. The second 
formula in the lemma follows from Stirling's formula. 

n 



4. A more precise asymptotic for the leading coefficients. 

In this section we assume that r = 0{k°'), with a < 2, unless otherwise stated. 

As a direct application of the Lagrange inversion formula, one can obtain a more 
precise expression for the saddle-point Ur^k over restricted ranges. This will allow 
us to write out series for u, U , Pk{u), and /i,„ which will give an explicit expansion 
for the saddle point approximation. 

Recall u is the unique positive number satisfying 






2k 



U+j 



y: ?^^=fc^ 



j=k+i 



U+J 



(4.1) 
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Since r < k", we know by lemma [4] that u > 2k for k large enough. So, we may 



apply Taylor expansions to the left side of equation (4.1 ). We therefore need to solve 

for u in the equation: X]m=i(~-*-)"'^^ ^ "> ^ ^ (^^"^ m = term cancels the fc^ on 
the other side of the equation). Substituting y = k/u the last equation becomes 

E™=i(-l)'"+'i^2/" = w- Now, trivially, 

0<p„,(A:)<(2fc)"+2, (4.2) 

so that the sum on the l.h.s. above is analytic in |y| < 1/2. Furthermore, it has 
non- vanishing derivative at y = 0, because pi{k) — k^. Using the bound on pm(fc), 
for X > sufficiently small, say < 1/100, one can find a unique y > 0, in say 
(0,1/10), such that 

Ec-ir+^^^z-^- (4-3) 

m— 1 

We can write down the Taylor series for y in terms of x by Lagrange inversion: 

y-^A„,(fc)a:™. (4.4) 

m— 1 

There are two formulas that are useful for finding the A's. We can either substi- 



tute the above formula for y into (4.3) and compare coefRcients of powers of x. 



Alternatively, we can use the well known formula (see [W], section 5.1): 

1 d' 
m\ dy 



^^^^) = -^^,:,^i^>^^yy 



(4.5) 



where Wfc(z) is given by Wife (z) := ^^^i(-l)™+i|^z"-i. Letting a„ = (-l)'"+V.(A:)/fc'"+2, 
we have 

Ai(fc) — 1 , A2(fc) — —02 , A3(fc) — 2al — a^ , A4(fc) — 50203 — 04 — 5al 
A5(fc) = 14 02** — 21 02^03 + 602O4 + 3o3^ — 05 

AeCfc) = — 42 02^ + 84 02*^03 — 2802^04 — 280203^ + 70205 + 70304 — og 
Xj{k) = 132 02^- 330 a2'*03 + 120 02^04 + 180 02^03^-36 02^05-72020304 
-1203^ + 80206 + 80305+404^-07 (4.6) 

Since y — k/u, and x — r/k"^, we get: 

l_r( iP-lr 22fc''-5fc2 + l / r \2 

M " fc3 l^^ ^ 6fc2 fc^ + 18P VfcsJ (-^ 

1357 k^ - 435 fc^ + 183 fc^ - 25 / r \ 3 4142 k^ - 1661 k*^ + 1083 fc"* - 359 fc^ + 35 



Vfc2j + 3240fc8 VpJ 



1080/c6 
58691 fcio - 28609 k^ + 27146 k^ - 14906 /c^ + 3283 fc^ - 245 / r \ 5 

45360fcio Vfc2 J 

888146 fci2 _ 506685 fc^o + 640353 k^ - 512890 fc^ + 201576 k^ - 32025 k^ + 1925 / r \ 6 



680400fci2 Vfc2 
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Because (4.4) is analytic in some neighbourhood of a; = 0, we have that there 
exists an 77 > such that the above converges geometrically for all r < rjk^, and, 
even though the coefRcients depend on k, does so uniformly in k. The latter is 



explained by the fact that the l.h.s. of (4.3 ) has, by (4.2 ), non-zero derivative in some 



neighbourhood, independent of k, oi y — 0, and that one can cover, independent 
of fc, some sufficiently small neighbourhood of x = by taking y small enough. 



Therefore, the inverse function specified in (4.4) is analytic about a; = in some 
neighbourhood that is independent of k. 

We can also write a series for u directly. On reciprocating; u — 

Using Maple, we compute: 



k 



i:™=iA„(fe)(^)'"- 



\m=l / 

- (A4 - 2 A2A3 + A2') x^ + (A2* - 3 As'As + 2 A2A4 + A3' - As) x^ 

- (A2^ - 2 A3A4 + 3 A2A3' - 4 A2^A3 + 3 A2'A4 - 2 A2A5 + Ag) a;^ 
+ (A2^ - 5 A2^A3 + 4 A2^A4 + 6 \2^\^^ - 3 Aa'As - 6 A2A3A4 - M^ 
+2A2A6 + 2A3A5 + A4'-A7)x6-... (4.8) 



Putting the above together gives the first few terms of u: 



_ k^ 7fc2-lr 5fc4-4fc2-l / r \2 

""~^ 6k^¥^ 36fc4 Vfc^J ^^■^> 



4 fc6 + 15 fc4 - 24 fc2 + 5 / r \ 3 59 /c^ - 152 fc6 - 114 A:4 + 232 k^ - 25 



U2J 



540fc6 V/cV 6480fc8 

92 fcio + 329 k^ - 2128 fc^ + 2302 k^ - 6Uk^ + 49 / r \5 
27216A?io Vpy 

3101 k^^ - 10620 fcio - 43827 k^ + 157640 k^ - 125649 fc'' + 20580 P - 1225 

1360800fci2 



A similar series for U can be derived from that of u. For u > 2fc — 1 we have: 



U 



(1 + j/u)^ "*" u .^ Jl 



'" J(2fc-j) 



^(-l)'"-'™p™(fc)7.-". (4.10) 
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Substituting u^^ ~ k^^ J2m=i ^m{k){r/k'^)™, the above equals 



oo 

r 



7k^ -1 r 5fc''-4fc2- 1 / r \2 
6fc2 F ^ 18fc4 Vfc2j 

151 fcS - 255 fc4 + 129 fc2 - 25 / r \ 3 187 fc*^ - 706 k^ + 168 fc^ + 386 fc^ - 35 / r \ 4 



V-^k^k^^ mi {¥) ^'-''^ 



VpJ ^ 3240fc8 VfcsJ 



1080fc6 
1373 /fci° - 9415 fcS + 18956 fc^ - 14540 fc^ + 3871 A:^ _ 245 

45360fcio V fc2 

7777 k^^ - 114000 fc^o + 158361 k^ + 190960 k^ - 279813 fc^ + 38640 k^ - 1925 
"^ 680400fci2 VP 

Next we work out the series for log(P/c(u)/(co(fc)u^ )). Divide the second formula 
in (1.31 by the leading coefficient Co(fc), pull out u*^ , expand each log(l + j/u) = 
j/u — j'^ /2u + . . ., and sum over j to get: 

log{P,{u)/{co{k)u'^')) = f: ^"'^""'^"^^^ . (4.12) 

?n— 1 

Substituting ^ = j: J2m=i ^m{k) ( j2 )™ the above equals 

m=l ^ ^ 

^ 7A:2-lr , 13fc4 - 2fc2 + l / r \2 583^*^ - 165fc'' + 147P - 25 / r y 



r 1 + 



(f)- 



12A:2 p 3gfc4 Vfc2y 2160fc6 \k^ ) 

1379/c«-428A;6 + 642fc'*-332fc2 + 35 / r \4 

6480fci (fc^) (^-^'^ 

3193 /fci° - 1393 k^ + 3178 k^ - 2542 fc^ + 637 fc^ - 49 / r \ 5 

18144A:io VpJ 

203849 fci2 - 100470 k^° + 307587 /c^ - 379420 k^ + 192639 fc-* - 31710 fc^ + 1925 / r \ 6 \ 



1360800fci2 

Finally, we develop series for ^m. These are deffiied by 1 + X]m=3A'm^™ == 
exp (X]^3 7n(w)^")! where 7n(M) are, from ( 2.27[ ), the Taylor coefficients of /(-ue*^). 



Vfc2J 



Substituting ue'" for z in /(z) = log(Pfc(z)) — (/c^ — r) log(z), we get, after Taylor 
expanding as in (4.121 and on expanding e*^ in its Taylor series, that 



m— 1 



n! 

771—1 



Article submitted to Royal Society 



18 



G.A. Hiary and M.O. Rubinstein^ 



Now each /Zm is given as a sum of products of the 7„(?i)'s, the specific expression 
obtained by exponentiating Y^^ ln{u)9^^: 



A^Tl 



E 

ii>0 



73 (m)^-' 74 (m)-'" 75 {uy'' 



i3!J4!J5! ■ 



(4.15) 



For example, /is = 78 + 7375 + 74/2. Notice that 7n(u) is r times a power series in 
r/fc^ with coefficients that are polynomials in 1/fc^. The latter follows from \j{k) 
being polynomial a polynomial in &„ = (— l)™+^pm(fc)/fc'"+^, and Pm{k) being 
a polynomial in k of degree m + 2 satisfying, as described in the introduction, 
Pm{k) = (— l)™p„j(— fc), so that the powers of k that appear 'go down by twos'. 
Therefore, each term in the above sum is of the form r° times a power series in 
r/fc^, with coefficients polynomial in 1/fc^, and a G Z, a < [w./3j . 

Substituting (IXl) (|4.11|) (|4.13|) (|4.15|) and (|4.14| into the saddle point for- 



mula (2.33), and taking the logarithm of the various series that appear as factors. 



suggests the following formula for Cr{k): 



Cr{k) 
Co(fe) 



U3r 



exp 



7r(r-l) r(r- l)(26r- 15) r(r - l)(583r2 - 715r + 183) 



12fc2 



144fc4 



6480fc6 



r{r - l)(2758r3 - 4499r2 + 463r + 1491) 
51840fc8 



•JjC^ 



fc2J 



(4.16) 



The j-ih term appearing in the exp is of the form qj(r)/k'^^ where qj{r) is a poly- 
nomial in r with rational coefficients and of degree j + 1. Furthermore, when r = 
or 1, Cr{k)/c^{k) = k^^/r\, and thus each term appearing in the exp must have 
numerator divisible by r{r — 1) so as to vanish at those two values of r. 

Below we develop our formula further, explaining with proof how it arises. Sub- 
stituting the various series as described does not immediately result in the above 
formula, as there are extraneous terms of the form am^„/r™fc2" along with various 
O terms that need to be considered. One also needs to identify the terms of Stir- 
ling's asymptotic formula as appearing from these substitutions so as to obtain the 



factor of r! in (4.16) 



These extraneous terms can be eliminated by comparing these terms with the 



terms that appear in equation (1.23), and considering different values of r. 

To begin, we have replaced the factors and terms with 'just an r and no k 
with a l/r!. The Pkiu)/{co{k)u''" 
of e"^ /r^ . The l/\/U contributes, from (|4.11|) a 1 



) contributes, from (4.13) and (4.9), a factor 



r, while the sum involving /i2m 
in ( 2.33[ ) contributes a factor whose initial terms are computed to be exp(— l/(12r)-|- 
l/(360r'^) + l/(1260r^) — ...). Together with the l/\/27r, we recognize these as 
Stirling's approximation to l/r!. By taking M in ( |2.33 ) sufficiently large, letting 
k ^ 00 and r — k'^ with e sufficiently small, and comparing with the asymptotic 
Cr{k)/cQ{k) — fc^''/r!(l + 0(fc2('^~i))) of Theorem[l] we have that these terms must 
coincide with Stirling's asymptotic series for l/r!. 

The extraneous terms can be eliminated by comparing with the full polynomial 
expansion of Cr(fc)/co (A;). From equations ( |1.20[ ) and (1.22), 

"-""-^^^ ' ' E 9Ar>' (4-17) 

l<i<(3r-2)/2 



k^^'coik) 



= 1 
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where z — 1/fc^, and gj{r) is a polynomial in r of degree < 2j (the bound on the 
degree of gj{r) follows from the inductive procedure for determining them). Denote 
the above polynomial in z by Fr{z). By taking its logarithm, we can write it in the 
form 

rlcr{k) 

k^'^coik) 



exp {qi{r)z + q2{r)z^ + qj{r)z^ + . . .) 



(4.18) 



where qj{r) is a polynomial in r. One can easily pass from the polynomial in (4.17) 



to an expression of the form (4.181, because the polynomial (4.17) does not vanish 
at z = 0, hence its logarithm is analytic there and may be expanded in a Taylor 



series. The constant term 1 in (4.17) implies that the Taylor series that appears in 
the exponent of ( 4.18[ ) has constant term 0, and thus begins with the qi{r)z term. 
Expanding out log Fr{z) using the Taylor series for log(l + w) — w~w'^/2 + w^ /3 — 



. . ., and comparing coefficients with those in the exponent of (4.18) shows that qj(r) 
is a polynomial in r of degree < 2j. Our goals are to show that qj{r) actually has 
degree j + 1 rather than < 2j, and to obtain a truncation bound for (4.18). 



To do so, we consider more carefully the result of substituting the various series 
into the saddle point formula, pulling out the various factors that comprise Stirling's 
asymptotic formula for r!, and moving everything else to the exp by taking the 
logarithm of the various factors: 

First, the series for log(Pfe(u)/(co(fc)u'^ )) is of the form r times a power series 
in r/fc^, whose jth coefficient is polynomial in 1/fc^ of degree j. 

Next, we consider the factor u^ . Pulling out, in ( |4.9[ ), fc^/r and moving the rest 
into the exp, we have that r log(u/(fc'^/r)) is of the same form as log{Pk{u)/u'' ) 
described in the previous paragraph. 

The series for 1/v C/ contributes, on pulling out a l/xA', and moving it to the 
exp as — log(f//r)/2, a power series in r/fc^ with coefficients in k as above. 

Finally, by the paragraph surrounding (4.15), fi2m/U"^ is equal to a sum whose 
terms are products of the 7„(u)'s, and hence expressible as power series in r/k^ 
with coefficients in k as above, with each term multiplied by a factor of the form 
1/r'^ with a > m — [2m/3\ . 
using, from (4.11), that U ^ 



Therefore, for M sufficiently large in (2.331 and also 
o(fc^), we get that the logarithm of the 



when r 



bracketed factor in (2.33) is, for any A > 0, of the form: 



-. oo 



(4.19) 



l<a<yl 



3=0 



where (Jj^aik) is a polynomial in 1/fc^ of degree j. 

Next, truncate all substituted power series in x = r/fc^ after, say, J = JA,a 
terms so that the remainder gets absorbed into the 0{r^^). This can be achieved 
because, by the paragraph following (4.7 1, all these power series can be shown 
to be analytic about a; = with their coefficients, expressed as polynomials in 
1/fc^, uniformly bounded in k. Hence each has some radius of convergence, and 
the cost of truncating each power series at J terms is 0{{ri'r/k^Y^^) for some 
rj' > depending on A. Finally, recall that r = 0{k°') with a < 2. Thus, taking J 
large enough (depending on A and a), we can ensure that the cost of truncating is 
Oir-"^). 
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Thus, collecting terms according to the power of 1/fc^, and pulling out the terms 
corresponding to Stirling's formula for r! as explained, we have obtained: 



.{k) 



k^''co{k) 



exp 






0{r 



-A\ 



(4.20) 



where Qj{A;r) is of the form 



~A+j<m<j + l 



(4.21) 



The largest power of r that appears is j ' + 1 rather than j because the power series 
for \og{Pk{u)/u'' ) and r log(M/(fc'^/r)) are multiplied by an extra r. The smallest 
power is > —A + j rather than —A because r™ /k'^^ < r™"^ can be asborbed into 
the 0{r~^) if m < -A + j. 



Formula ( |4.20 l essentially gives, as a function of z = l/k^, the Taylor expansion 
of \og{rlcr{k)/{k^^ca{k))), and comparing with (4.18), we would like to conclude 
that qj{r) = Qj{A]r), from which we would have that qj{r) is a polynomial of 
degree j + 1, rather than of degree 2j. However, the presence of extraneous terms 



and the remainder term in (4.20) necessitates a more careful comparison, and we 



will only conclude that the polynomial part of both coincide for given j and all A 
sufficiently large. 



From equation (2.14), we have that r\cr{k) / {k^^ CQ{k)) is asymptotically 1 when 
r = O^k"), for a < 1. Therefore its logarithm is analytic in z = 1/k'^ in that region. 



and its Taylor series in (4.18) converges in at least that region. 



For r — 0{k°'), a < 1 we can truncate (4.18 1 after sufficiently many terms so as 
to also have remainder 0{r^ ). By adjusting our value of J upwards, we can assume 
that both series in (4.18) and (|4.20l truncate at J. We thus have, on comparing 



them: J2j=i 



J 9,(r)-QjT^?FF 



0{r- 



^ j^ -p3 — <^v' )■ We can use this to show that, for any given j, 

and all sufficiently large A, that qj{r) — Qj{A; r). We have established that qj{r) is 



a polynomial of degree < 2j in r, and that Qj{A; r) has the form (4.21 ). Therefore, 
the above can be written as: X]i=i ^-A+j<m.<2j ^"ij" W = 0{r^^). For each choice 
of A, the l.h.s. above consists of finitely many terms. If we let r ^ k'^ , with a < 1 
and a irrational so as to guarantee that all the exponents ma — 2j are distinct, 
then each term above must individually satisfy the bound: dm.jk"^°'~^^ = 0{k~^"). 



Therefore, for fixed j, m > and all large enough A, we have that d„ 



0, 



qj(r) coincides with the polynomial portion of Qj{A; r). This shows that qj{r) has 
degree j + 1- 

Next we obtain a bound on the rate of growth of qj{r). We have seen qj{r) 
coincides with the polynomial portion of Qj{A; r) if we take sufficiently many terms 
m in (2.33), and have observed that fi-zm/U™' has a factor of l/r" with a = m — 
[2m/3\ > m/3. Thus, given j, in order for the various substituted power series to 
capture the degree j + 1 polynomial (/^(r), we need to take, in (2.33), |[(M— 1)/2J > 
j + 1, so M > 6(j + 1) suffices. Thus we let M = 6{j + 1) , and consider: 



log 



Pk{u) 



fc3'-co(fc)Mfe' y/2nU 



'1^'^ 2'"r(m + 1/2) ^2. 

m=2 



/ttU''' 



(4.22) 
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As before, we substitute the power series in r/fc^ for the various factors and terms 
of this expression, and letting z = l/k^, express it as 

oo 

5^Q,,,(r)z\ (4.23) 

2=1 

We will obtain a bound for Qij{r) and then deduce a bound on qj{r) by bounding 
the polynomial part of Qj,j{r). To do so, we regard the above expansion as a 
function of a complex variable r and z and apply Cauchy's estimate twice- first to 
bound Qij{r), and then to bound its polynomial part. 

From the power series expansions for log(Pfe(M)/(co(fc)u*^ )), r\og{u/{k^/r)), 
and — log(C//r)/2, we have that: there exists a > such that for all \z\ < a/\r\ 
(think of this as r/fc^ < a), such that 

log (^^|i^) ^0(0. (4.24) 

Furthermore, in the same kind of region, \z\ < a/|r|, we have f/ — r2(|r|), and 
[i-m = Odrl™'"^), uniformly in m. Combining this with the elementary estimate, 
r(m + 1/2) < TO™V^, gives 

r(^ + l/2),.™,_^(^.^^g^.^^ (4.25) 




f^W' 



so long as \r\ is bounded away from 0. Therefore, applying Cauchy's estimate 
to (4.23[), taking the circle \z\ = a/|r| as the contour of integration, gives Qi.j{r) = 



0{j\og{j)r{r/ay). In particular, Qj,j{r) = 0{j\og{j)r{r/ay). 

To bound the polynomial part qj{r) of Qjj{r) we take the integral 



j+i 



•^^ rn=l 

with C the circle \w\ = 2. Applying the above estimate for Qjj{r) gives qj{r) = 
0{jlog{j)r{rXy), for some A > 0. 



Presumably, there should be a more direct and combinatorial derivation of (4.16 ) 



One can verify that expanding (4.161 using the Taylor series for exp does produce 



the first few terms of (1.231. In fact, by comparing these two expressions, one can 



more easily calculate the polynomials qj{r) from the terms in (1.23) 



5. A maximal co efficient 

In this section we obtain the estimate for the maximal Cr{k) given in Theorem pi 

The Cr(fc)'s are the coefficients of a polynomial with negative roots. This im- 
plies the sequence {cr{k)} is unimodal (see jW]). That is, the sequence {cr{k)} 
monotonically increases until it reaches a maximum, then it monotonically de- 
creases. Or there exists s € Z"*" such that co(fc) < ci{k) < ■■■ < Cs{k), and 
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Cs{k) > Cs^i{k) > ■ ■ • > cj.2(fc). The idea is to combine unimodality with a saddle- 
point estimate for Cr+i{k) — Cr{k). Adapting the proof of lemma pi we consider 

A. := c.+i(fc) -cr{k)^^ f ^' ',PS^i'^ dz , (5.1) 

where C :— {ue*^ : — tt < 6 < tt}. Choose u :— Uj.,k to be the saddle-point corre- 
sponding to Cr{k); that is, u is the positive number satisfying 

uY^+uY'^^^^^e-r. (5.2) 

3 = 1 3=k+l 

Let A be defined as in ( |l.8| ), so A = log4-|-l/(2/c) + 0(l/fc^). An examination of 
the results in Section 2 reveals a maximal coefficient Cr{k) occurs when the saddle 
point Ur^k ~ 1- By lemma Wl this happens when k/^ — r ^ kA =^ r « fc^ — /x. 
Guided by this, our plan is to show, for some p > and all k sufficiently large, 
that there exists integers ri e (fc^ — ^ + p\og{k)'^/k, fc^ — ^ + 1 + plog(fc)^/fc] and 
r2 G [fc^ — ^—l — plog(/i;)^/fc,fc^ — ^ — plog(fc)^/A:), such that A^i is strictly negative, 
and Ar2 is strictly positive. Combined with the unimodality of the sequence {cr{k)}, 
this shows Cr{k) is maximal for some r e [k^ — fi— plog{k)'^ / k, k'^—fi+l+p\og{k)'^ /k], 
and is not maximal outside of that interval. To this end, consider the behaviour of 
Ar for 

re[k^-kA^2 log{kf , fc^ - fcA + 2 log(fc)2] . (5.3) 

Under assumption ( |5.3| , it follows by lemma l4J 

k^ — r /log(fc)^ 



UrM = -^-r- [1 + 0(log(fc)/A:)] = 1 + 



kA ' ^ ^' " " \ k 

U := Ur,k^{k^-r)[l + O{\0g{k)/k)]^kA + 0i\og{k)''). (5.4) 

Using the same procedure as in the proof of lemma Isl and choosing 6 — 1/1000, 



M = 13, in expression (2.301, one obtains 



where e = t/i/iooo-i/2^ ^^^^ ^^ ^ Oj{U^^/^i). Replacing ue*^ with u + iue~ue^/2 + 
0{e^) yields 



/ {u + iue-ue^/2-i)\i + Y ^J']^' \ 



^'^" rf^ + O 77^ .(5.6) 



The integral vanishes for odd powers of 9. Also, the contribution of some the terms 
can be absorbed into the remainder. For example 

/ {u-l)\l + Y^l,ei\e''i'"de^ I (u-i)e-*^'d0 + o(c/-^/2) .(5.7) 
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where we made use of u — 1 — 0((log k)/k), which holds according to (5.4), the 



bound on fj,j, and also that odd powers of 6 integrate to 0. After such reductions, 
the polynomial that survives is m — 1 — ufP' jl + iu^^d'^. So, when the domain of 
integration is extended to all of M, we arrive at 



Define 



~2^ 



{u-i~ ue^/2 + iufi3e'^)e-i''^ de + o 



\u^ 



/2Tr J- 



n TT—ni — 1/2 I -- „ 



/27r 



(5.8) 



.90 



C/l/2' 



51 



[/3/2' 



91 



1/5/2 



Then, 



/27rA 



-^ = {u - l)go - M.gi/2 + im3.92 + O (u~^^^^ 



(5.9) 



(5.10) 



Thus, for k greater than some absolute constant and for r in the range (5.3) 

u-l + R, 



V^ttU A^ 



=/(") 



(5.11) 



where 



R 



u 3iuj:i 



+ 0{U-^). 



(5.12) 



2U C/2 

Next we obtain a more precise estimate for /13 = 73. Differentiating f{ue^^) 
three times with respect to 9 and setting 6 — gives: 



73 



lU 

31 



E 

i=i 
iu 
3[ 



f{u-j) 



2k 



E 

j=k+i 



{2k-j)j{u-j) 



(fclog4-4Mlogfc + 0(l)). 



(5.13) 



Now, from (|5.3[) and (1.8), we get: 73 



iuU 



0(log(/c)^). Therefore, plugging 
into (5.12) and simplifying gives R = -^{u— \)^-0{U ~'^ {log kY) (we have dropped 
the 0(1/C/^) which is subsumed by the 0(C/"^(logfc)2)). Hence: 



^^^§^ = {u-m + ^)+o{u~-iiogkf). 



(5.14) 



The sign of A^ is, therefore, essentially determined by that of u — 1. To obtain the 



maximal Cr(fc) we develop, below, a stronger estimate for u than (5.4 1. 
Recall u is the positive number satisfying 



h{u) 



k" 



where h{x) :— x \^ - 



] 



j=i 



] +x 



y. 2k-j 



j=fc+i 



J +x 



(5.15) 
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Also, one may easily verify U — uh'{u). Expanding in a Taylor series, 
h{l + x)^ h{l) + h'{\)x + 0{h"{l)x^) , 



(5.16) 



for \x\ < 1/2 say. And by direct calculation, h{l) 



Ek j \r^2k 2k— j 

3=1 7+T ^ l^j=k+i 7+r 

^2k {2k-])j _ 

^ 2^3=k+l (j + l)2 



= A; log 4 



kA - \og{kl2) - 7 + 0{l/k), h'{l) = E;=i (liTF 
O(logfc), and h"{l) ^ O(logfc). 

Let p > 0, and choose d G (plog(A:)^/fc, 1 + plog(A:)^/fc] so that h{l) — d 
is an integer. Consider the sadd le po int u^ corresponding to C}^2_i^(^i-^j^.^{k)] i.e. 
h{u^) = h{l) — d. By expansion (5.161, the point u^ satisfies u~ = 1 — d/h'{l) + 
0((l-u^)2logfc//i'(l)), Since fc-^ - h{l) + d falls in the range (O, it follows 



from estimate (5.4) that 1 — u^ = 0(log k/k). Consequently, u^ = 1 — d/h'[l) + 
0(log(fc)'^/fc'^). Furthermore, by (5.4), the corresponding U satisfies U ^ kA. Thus (5.14) 
becomes 

"^ 0{k-^{\ogkf). (5.17) 



\/27rC/Afe2_ 



fe2-h(l)+d 



=/(«) 



;.'(i) 



So, choosing p sufficiently large, we have, for all sufficiently large fc, that lS.}^2_f^(j-^^^ 
is strictly negative. Similarly, if one chooses d € [—1 — plog(fc)^/A:, — plog(fc)^/fc) 
we obtain a strict inequality in the opposite direction. Theorem [2] follows. 

6. Numerical verifications 

TableII|compares different approximations to Cr(7). Recall that 5^(7) := Cr(7)/co(7). 
Column 3 illustrates that, for smaller r, the asymptotics of hr{k) are simply 

described by fc''(^), while column 7, which compares Cr(fc) to (^,2_^)(A/A:)'' ~'', 
works best for r near k^ . 

Column 4 shows the ratio of br{k) to the explicit formula given in (4.16), trun- 
cating the terms in the exp at those that are depicted, namely terms up to Ijk^ . 

Column 5 compares the ratio of Cr(7) to the saddle point formula (2.18), without 
the O term. It approximates well throughout, except near the two tails. 

Finally, column 6 shows the performance of the uniform approximation, depict- 
ing Cr(7) divided by the r.h.s. of ( |1.5| , without the O-term. The uniform asymptotic 
performs surprisingly well across all ranges, and the ratio seems to equal l+0(l/fc^) 
uniformly in r, rather than 1 -I- O (logfc/fc^ ^'^) as proven in Lemma 5 



This is somewhat surprising because (1.5 1 was designed to interpolate three 
formulas co{k)k^'^ /r\, (2.18), and {kAy /rl, across three asymptotic ranges. Yet, for 
every r, it seems to approximate better than any of these. Presumably the uniform 
approximation correctly captures lower order terms especially in the tails where the 
improvement is substantial. 

This is also illustrated by Figure 1 which compares the performance of the 
uniform asymptotic to the saddle point asymptotic, for k = 100. We see the uniform 
asymptotic agreeing to within 1/fc^ across all r. 

References 

[CFKRSl] J.B. Conrey, D.W. Farmer, J. P. Keating, M.O. Rubinstein, N.C. Snaith, Integral 
moments of L -functions, Proc. London Math. Soc. (3) 91 (2005), no. 1, 33-104. 



Article submitted to Royal Society 



Uniform asymptotics of the coefficients of unitary moment polynomials 25 



[CFKRS2] J.B. Conrey, D.W. Farmer, J.P. Keating, M.O. Rubinstein, N.C. Snaith, Lower 

order terms in the full moment conjecture for the Riemann zeta function, Journal of 

Number Theory, Volume 128, Issue 6, June 2008, 1516-1554. 
[CGI] J.B. Conrey, A. Ghosh, Mean values of the Riemann zeta function, Mathematika 31 

(1984) 159-161. 
[CG2] J.B. Conrey, A. Ghosh, A conjecture for the sixth power moment of the Riemann 

zeta function, Int. Math. Res. Not., 15 (1998), 775-780. 
[CGo] J.B. Conrey, S.M. Gonek, High Moments of the Riemann Zeta- function, Duke Math. 

J., V. 107, No. 3 (2001), 577-604. 
[CRT] R. Chelluri, L.B. Richmond, N.M. Temme, Asymptotic estimates for generalized 

Stirling numbers, Analysis (Munich) 20 (2000), no. 1, 1-13. 
[DB] N.G. de Bruijn, Asymptotic methods in analysis, Dover Publications, Inc., New York, 

1981. 
[FGH] D. Farmer, S. Gonek, C. Hughes, The maximum size of L-functions, Journal fiir die 

Reine und Ange. Math., 609 (2007), 215-236. 
[HR] G. Hiary, M. Rubinstein Uniform asymptotics of the coefficients of the moment poly- 
nomials for the Riemann zeta function, preprint. 
[KS] J.P. Keating, N.C. Snaith, Random matrix theory and (^(1/2 + ii), Comm. Math. Phys. 

214 (2000), no. 1, 57-89. 
[MW] L. Moser, M. Wyman, Asymptotic development of the Stirling numbers of the first 

kind, J. London Math. Soc. 33 1958 133-146. 
[W] H.S. Will, generatingfunctionology. Third edition. A K Peters, Ltd., Wellesley, MA, 

2006. x-l-245 pp. 



Article submitted to Royal Society 



26 



G.A. Hiary and M.O. Rubinstein^ 



r 


br{7) 


6,(7) 


1H6| 


1H8} 
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0.92753 


46 


1.51291e+41 


0.0109627 


0.01840875972 


0.97341 


1.00254 


0.962089 


47 


5.26306e+40 


0.00853539 


0.00698803733 


0.960117 


1.00242 


0.986828 


48 


1.18076e+40 


0.0065654 


0.001898534367 


0.922532 


1.00217 


1 


49 


1.28039e+39 


0.00498356 


0.0002772583691 


2.05827e-05 


Inf 


1 



Table 1. A comparison of five formulas for the asymptotics o/c,.(7), as explained in 

Section\^ 
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k 


r 


k^ -,1 


r — k'^ + 11 


k 


r 


k^ ~ fi 


r - k^ + fi 


2 


3 


2.5833333 


0.41666667 


51 


2535 


2534.6271 


0.37291935 


3 


7 


6.5166667 


0.48333333 


52 


2637 


2636.26 


0.74002916 


4 


13 


12.372619 


0.62738095 


53 


2740 


2739.8925 


0.1075001 


5 


20 


20.181349 


-0.18134921 


54 


2846 


2845.5247 


0.47531881 


6 


30 


29.958261 


0.041738817 


55 


2954 


2953.1565 


0.84347267 


7 


42 


41.712279 


0.28772061 


56 


3063 


3062.7881 


0.21194975 


8 


56 


55.449036 


0.55096431 


57 


3175 


3174.4193 


0.58073872 


9 


72 


71.172312 


0.82768841 


58 


3289 


3288.0502 


0.94982886 


10 


89 


88.884769 


0.11523121 


59 


3404 


3403.6808 


0.31920997 


11 


109 


108.58835 


0.41164848 


60 


3522 


3521.3111 


0.68887239 


12 


131 


130.28452 


0.71547681 


61 


3641 


3640.9412 


0.0588069 


13 


154 


153.97441 


0.025587194 


62 


3763 


3762.571 


0.42900475 


14 


180 


179.65891 


0.34109033 


63 


3887 


3886.2005 


0.79945758 


15 


208 


207.33873 


0.66127328 


64 


4012 


4011.8298 


0.17015746 


16 


237 


237.01444 


-0.01444432 


65 


4140 


4139.4589 


0.54109679 


17 


269 


268.68654 


0.31345885 


66 


4270 


4269.0877 


0.91226833 


18 


303 


302.35542 


0.64458331 


67 


4401 


4400.7163 


0.28366517 


19 


338 


338.02141 


-0.021407767 


68 


4535 


4534.3447 


0.65528068 


20 


376 


375.6848 


0.315199 


69 


4670 


4669.9729 


0.02710854 


21 


416 


415.34584 


0.65415768 


70 


4808 


4807.6009 


0.3991427 


22 


457 


457.00474 


-0.0047443102 


71 


4948 


4947.2286 


0.77137735 


23 


501 


500.66169 


0.33830803 


72 


5089 


5088.8562 


0.14380693 


24 


547 


546.31685 


0.68315269 


73 


5233 


5232.4836 


0.51642612 


25 


594 


593.97035 


0.029647012 


74 


5379 


5378.1108 


0.88922979 


26 


644 


643.62234 


0.37766471 


75 


5526 


5525.7378 


0.26221304 


27 


696 


695.27291 


0.72709347 


76 


5676 


5675.3646 


0.63537115 


28 


749 


748.92217 


0.077832963 


77 


5827 


5826.9913 


0.0086995789 


29 


805 


804.57021 


0.42979318 


78 


5981 


5980.6178 


0.38219397 


30 


863 


862.21711 


0.78289307 


79 


6137 


6136.2441 


0.75585013 


31 


922 


921.86294 


0.13705943 


80 


6294 


6293.8703 


0.12966403 


32 


984 


983.50777 


0.49222584 


81 


6454 


6453.4964 


0.50363176 


33 


1048 


1047.1517 


0.84833194 


82 


6616 


6615.1223 


0.87774959 


34 


1113 


1112.7947 


0.20532266 


83 


6779 


6778.748 


0.2520139 


35 


1181 


1180.4369 


0.56314762 


84 


6945 


6944.3736 


0.62642122 


36 


1251 


1250.0782 


0.92176065 


85 


7112 


7111.999 


0.00096817198 


37 


1322 


1321.7189 


0.28111928 


86 


7282 


7281.6243 


0.37565153 


38 


1396 


1395.3588 


0.64118441 


87 


7454 


7453.2495 


0.75046815 


39 


1471 


1470.9981 


0.0019199218 


88 


7627 


7626.8746 


0.12541502 


40 


1549 


1548.6367 


0.3632924 


89 


7803 


7802.4995 


0.5004892 


41 


1629 


1628.2747 


0.72527087 


90 


7981 


7980.1243 


0.87568788 


42 


1710 


1709.9122 


0.087826568 


91 


8160 


8159.749 


0.25100833 


43 


1794 


1793.5491 


0.45093271 


92 


8342 


8341.3736 


0.62644789 


44 


1880 


1879.1854 


0.81456437 


93 


8525 


8524.998 


0.0020040182 


45 


1967 


1966.8213 


0.17869823 


94 


8711 


8710.6223 


0.37767423 


46 


2057 


2056.4567 


0.54331253 


95 


8899 


8898.2465 


0.75345612 


47 


2149 


2148.0916 


0.90838688 


96 


9088 


9087.8707 


0.12934738 


48 


2242 


2241.7261 


0.27390216 


97 


9280 


9279.4947 


0.50534575 


49 


2338 


2337.3602 


0.63984041 


98 


9474 


9473.1186 


0.88144904 


50 


2435 


2434.9938 


0.006184772 


99 


9669 


9668.7423 


0.25765515 



Table 2. A comparison of the maximal Cr(fe) for 2 < k < 99, to the value k^ — /i of 
TheoremY^ The Theorem states that, there exists p > such that, for all k sufficiently 
large, their difference lies in [— plogfc^/fc, 1 + plogfc^/fc]. 
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Figure 1. We compare three approximations for Cr(lOO), < r < 10000. The long-dashed 
line depicts the ratio of Cr(lOO) to the r.h.s. of the uniform approximation ( |1.5[ ), while the 
solid line shows the ratio of Cr-(IOO) to the r.h.s. of the saddle point approximation (2.181, 
in both cases without the O-term. The short-dashed line plots the ratio of Cr(100)/co(100) 
to the formula shown in ( |4.16[ ), taking the first 4 terms in the exp, i.e. dropping terms 
with a l/fc^" or higher. 
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